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SUMMARY 

Finite-element solution technology developed for use in conjunction with 
advanced viscoplastic models is described. The development of such solution 
technology is necessary for performing stress/life analyses of engineering 
structural problems where the complex geometries and loadings make the con- 
ventional analytical solutions difficult. The versatility of the solution 
technology is demonstrated by applying it to viscoplastic models possessing 
different mathematical structures and encompassing isotropic and anisotropic 
materials. The computational results qualitatively replicate deformation 
behavior observed in experiments on prototypical structural components. 


ORIGIN CQKTMNS 


COLOR ^ s ' 


iTRKUOH^ 


INTRODUCTION 

The drive to enhance and improve the high-temperature performance of 
structural components used in aerospace and nuclear industries has created a 
need for the development of realistic viscoplastic models, accompanied by 
appropriate solution technologies for stress/life analyses of these components. 
The observed interaction between creep and plastic deformation at high tempera- 
tures has led to the development of a number of viscoplastic models. These 
models treat all Inelastic strain as a single time-dependent quantity, and 
thus automatically Include creep, relaxation, and plasticity Interactions. 

Viscoplastic models are intended to provide realistic descriptions of 
high-temperature, time-dependent, inelastic behavior of materials. Their 
mathematical structure is, however, very complex. The highly nonlinear and 
mathematically "stiff" nature of the constitutive equations makes closed-form 
analytical solutions virtually impossible. It is, therefore, of the utmost im- 
portance that suitable finite-element or other numerical-solution technologies 
be developed to make these models adaptable for realistic structural and life 
analyses of these components. 

Researchers at the NASA Lewis Research Center have undertaken this impor- 
tant and challenging task. As a result of concerted efforts at NASA Lewis dur- 
ing the last few years, several such solution technologies (in conjunction with 
the finite-element program MARC and other nonlinear structural analysis codes) 
have been developed and successfully applied to the solution of a number of 
uniaxial and multiaxial problems. 

This report illustrates the application of some of these solution technol- 
ogies with regard to the viscoplastic models developed by Robinson, Freed, and 
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Dame and Stouffer (refs. 1 to 3). These models have different mathematical 
structures and encompass both isotropic and anisotropic materials. The intent 
of this work is to familiarize the researchers and designers in industry with 
the applicability of these technologies to structural engineering problems. It 
is expected that the results from the present work will stimulate the use of 
realistic viscoplastic models for better and more rational component designs. 


FINITE-ELEMENT SOLUTION TECHNOLOGY 

The finite-element analyses for the uniaxial and multiaxial problems were 
performed by using the finite-element program MARC (ref. 4). The program has 
been developed expressly for nonlinear structural analyses and is well equipped 
with sophisticated computational algorithms and advanced finite-element formu- 
lations. The details of the finite-element solution technology developed for 
use in conjunction with Robinson's model are given in reference 5. The tech- 
nologies for the other models follow essentially the same structure and can 
easily be developed. The nonlinear constitutive relationships of the visco- 
plastic .model s are implemented through the user subroutine HYPELA. The stiff 
nature of the constitutive equations requires smart and efficient time- 
integration strategies for their integration. The integration strategy em- 
ployed for the present computations is described in the next section. 


INTEGRATION STRATEGY 

The nonlinear and mathematically stiff nature of the constitutive equa- 
tions of v.iscoplastic models calls for special strategies for their integra- 
tion. This requirement has long been recognized by researchers. As a result, 
numerous integration strategies that employ either the implicit integration 
methods (refs. 6 to 12) or the explicit integration methods (refs. 6, 10, 

A 13-15) have been proposed. Since the use of an explicit integration strategy 
does not require evaluating or inverting a Jacobian matrix, the Euler forward 
method with an automatic time-step size was used in the present work. A de- 
tailed account of various explicit automatic time-step integration strategies 
is given in references 13 and 15. 

The Integration strategy used herein can be summarized as follows: 

Consider the equation 

^ = F(y,t) (1) 

For the solution at time t + At to be obtained from the known solution at 
time t, the following steps are required. 

(1) Choose a suitable time step (such as At). Denote the value of y 
obtained with this time step as Yj (= yt+At)- 

(2) Halve the time step At. Denote the value of y (obtained in two 
steps) as Y n (= y t+ At>- 
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(3) Define the tolerance as 


Tol 


1*1 - y ni 

l y Ill 


The upper and lower limits for Tol are prescribed and denoted by Tol u and 
Tol q, , respectively. The strategy now runs as follows: 

(4) If Tol > Tol u , replace At by At/2. Go to step (1). Recompute 
Tol by repeating steps (1) to (3). Repeat the procedure until a time step 
At such that Tol < Tol u is obtained. Accept the corresponding value of 
y as the value yt+At- 

(5) Compare Tol with Tol g, . If ToIq, > Tol, double the time step as 
used in step (4). If Tolg, < Tol, retain the time step as used in step (4) 
for the next step calculations. 

(6) Continue the computations with the time step rendered by step (5), 
and so on. 

The strategy is easy to implement and works well even for problems with 
complex geometries and loadings. It can easily be generalized for a system of 
equations by using concepts (such as defining the error norms) similar to those 
given in references 13 and '15. 


ROBINSON'S MODEL 

The viscoplastic model developed by Robinson utilizes the concept of yield 
surface. The form of model used in the present work has only one state varia- 
ble, namely, the back stress variable. The yield stress variable is assumed 
to be constant. The flow law, evolutionary law, and material constants for 
2-l/4Cr-1Mo steel can be found in reference 1. The viscoplastic model was 
implemented in the MARC program and the details of implementation are available 
in reference 5. Several uniaxial and multiaxial problems were analyzed to dem- 
onstrate the feasibility and applicability of the implementation as a useful 
structural analysis tool. Some representative results are presented below. 

A comparison of MARC and experimental hysteresis loops for the alloy 
2- 1 /4Cr- 1 Mo steel at two strain rates is presented in figure 1. Excellent 
agreement is observed between the MARC and experimental loops. Such uniaxial 
calculations (both here and in the case of other viscoplastic models) and their 
comparison with the experimental data are done in order to validate the finite- 
element implementation of the model and to gain confidence in applying it to 
complex and computationally intensive problems. 

As an application of the finite-element solution technology developed at 
NASA Lewis to an engineering structural problem, a nonlinear structural analy- 
sis was performed for an actively cooled cowl lip (the name given to the struc- 
tural leading edge of an engine inlet of a proposed hypersonic aircraft). This 
structural component is subject to severe thermal loadings and gradients during 
simulated flight. The three-dimensional finite-element model (shown in fig. 2) 
and the calculated steady-state temperatures on the cowl lip were taken from 
the work of Melis and Gladden (ref. 16). The material of the cowl lip was 
assumed to be a copper-based alloy NARloy-Z. The constants for this material 
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were obtained from Robinson and Arnold (ref. 17). The thermal loading cycle 
used in numerical computations is shown in figure 3. 

Figures 4 to 6 show the stress and strain distributions in the cowl lip. 
The stress and total strain distributions along the edge of the cowl lip at 
0.75 sec are shown in figures 4 and 5, respectively. The total (mechanical) 
strain in figure 5 is the sum of elastic and inelastic strains. From figure 5 
the largest total strain is observed to occur along the edge of the cowl lip. 
Figure 6 displays the stresses at 2.25 sec. It is interesting to note that 
the elastic-viscoplastic analysis is capable of capturing the redistribution 
of stress during even a small period of steady-state thermal loading, that is, 
between 0.75 and 2.25 sec. This redistribution of stress with time is evident 
on comparing figure 4 with figure 6. The elastic-plastic-creep analysis in 
reference 18 was unable to depict this redistribution of stress with time. 


FREED'S VISCOPLASTIC MODEL 

Freed's model incorporates two internal state variables, namely, the back 
stress and drag stress. The tensorial back stress accounts for strain-induced 
kinematic hardening, whereas the isotropic hardening in the model is included 
via the scalar drag stress. A small displacement and a small strain formula- 
tion are employed. It is assumed that there is no coupling between the static 
and dynamic recovery terms in the evolutionary equations for the state varia- 
bles. This imparts a simpler mathematical framework for the model. The flow 
law, evolutionary laws, and values of the constants for copper for this model 
can be found in reference 2. The model was implemented in the MARC program, 
and several uniaxial problems and a cylindrical thrust chamber nozzle problem 
were analyzed by using this finite-element implementation. The uniaxial load- 
ings shown in figure 7 consisted of mechanical and thermal (both isothermal 
and noni sothermal ) loadings. Noni sothermal loadings included both in-phase 
and out-of-phase loadings. 

Stabilized hysteresis loops, calculated from the model, for isothermal and 
noni sothermal (both in-phase and out-of-phase) loadings are plotted in fig- 
ures 8 to 10. The strain range is 2.0 percent and the strain rate is 0.001/sec 
for isothermal loops. For noni sothermal loops, the mechanical strain range is 
0.8 percent and the rate is 1 ,5xl0 - ^/sec. To facilitate comparison, the exper- 
imental hysteresis loops are also plotted in these figures. A comparison of 
these figures reveals good agreement between the experimental and predicted 
stabilized hysteresis loops obtained with the finite-element implementation. 
This gives confidence to apply the implementation to a multiaxial nozzle 
probl em. 

The finite-element model of the cylindrical rocket nozzle shown in fig- 
ure 11 and the calculated temperatures and pressure loadings were taken from 
reference 19. Generalized plane-strain isoparametric elements were used to 
model the smallest repeating segment of the cylinder wall. Time-varying nodal 
temperatures and pressure loadings (see ref. 19) were applied to compute cham- 
ber wall deformation histories for a number of firing cycles. The nonlinear 
variations in temperature-dependent material properties were accounted for in 
the computations. 

The deformed geometries of the component after the first and sixth firing 
cycles are plotted in figure 12. A magnification factor of about 15 was 
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applied to the deformed shape to facilitate the visual interpretation of the 
results. Also shown Is the deformed shape of the component as observed in a 
cyclic experiment carried to 393 cycles. It is encouraging to note that the 
analysis qualitatively predicts the progressively deformed shape that result's 
in the "dog-house" effect observed In the experiments. 


STOUFFER'S SINGLE-CRYSTAL MODEL 

I 

Owing to the presence of grain boundaries, conventionally cast polycrys- 
talline superalloys are susceptible to transverse grain boundary oxidation, 
corrosion and creep deformation, and subsequent cracking. The absence of grain 
boundaries in single-crystal alloys such as Rene-N4, MAR-M 247, PWA 1480, and 
CMSX has made them useful in gas turbine engines. A time-dependent, micro- 
mechanics, crystallographic model was developed by Dame and Stouffer (ref. 3) 
to characterize the inelastic behavior of single-crystal alloys. The model 
contains only the isotropic state variable. To make the model applicable in 
the analysis and design of single-crystal components, it was implemented in 
the MARC program. Several uniaxial problems, including creep, relaxation, and 
tensile and cyclic loadings, were analyzed. 

Predicted and experimental stress-strain curves for the single-crystal 
alloy Rene-N4 at 1400 °F (760 °C) are plotted in figure 13. The crystal orien- 
tations are (1 0 0), <1 1 0) and (930). Good agreement between predictions 
using the MARC program and experimental results was observed. 

Figure 14 exhibits creep curves at 1400 °F (760 °C) in the (1 0 0) direc- 
tion. The curves are plotted for three different loadings of 621 MPa (90 ksi), 
655 MPa (95 ksi), and 758 MPa (110 ksi). The experimental curves are also 
plotted to facilitate comparison. Again, good agreement between the predicted 
and experimental results was obtained. 

A comparison of predicted relaxation curves from the present work and the 
work of Dame and Stouffer (ref. 3) is presented in figure 15. The crystal 
orientations are (1 0 0) and (1 10). The experimental relaxation data for 
Rene-N4 were unavailable. It is seen that the predicted curves from these two 
works compare very well. 

Figure 16 displays a comparison of MARC calculated and experimental hys- 
teresis loops in the (1 1 0) direction. The strain range is ±0.6 percent and 
the temperature is 1400 °F (760 °C) . Excellent agreement between the predicted 
and experimental hysteresis loops confirmed the correct finite-element imple- 
mentation of the single-crystal constitutive model. It gives the confidence 
to apply this finite-element implementation to more complex loading conditions 
of more complex geometries. The results of such analyses will be reported 
subsequently. 


CONCLUSIONS 

Numerical-solution technology employing the finite-element method was 
developed to perform the nonlinear structural analyses with viscoplastic 
models. The technology is Illustrated for three models but is general in 
nature. Several uniaxial problems and a mu 1 1 1 ax 1 al problem were analyzed. 
Good agreement between the experimental and predicted hysteresis loops for 
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uniaxial thermomechanical loadings confirms the correct finite-element imple- 
mentation of the viscoplastic models. Finite-element analyses were performed 
for the multiaxial problems by using the "validated" implementation. It is 
seen that viscoplastic models are capable of describing thermal /mechani cal 
phenomena observed in the experiments. The versatility of the finite-element 
solution technology developed herein is demonstrated by successful application 
to three different types of viscoplastic models. 
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FIGURE 1. - STABILIZED HYSTERESIS LOOPS FOR 2-1/M Cr-1 Mo 
STEEL AT 538 °C (1000 °F). EXPERIMENTAL RESULTS TAKEN 
FROM ROBINSON AND SWINDEMAN (REF. 1). 


FIGURE 2. - COWL LIP FINITE -ELEMENT MODEL. ELEMENTS, 329M; 
NODES, M760. 



FIGURE 3. - THERMAL LOADING CYCLE FOR COWL LIP 
PROBLEM. 
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STRESS 


ksi 16.9 11.1 5.33 -.466 -6.26 -12.1 -17.9 -23.7 



116.5 76.5 36.8 -3.21 -432 -63.4 -123.4 -163.4 


MPa 



2 X 


FIGURE 4. - STRESS IN Z-DIRECTION IN COWL LIP; ELASTIC-VISCOPLASTIC ANALYSIS (ROBINSON’S MODEL). TIME. 0.75 sec; CROSS FLOW, 
GH 2 COOLANT. 
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FIGURE 5. - TOTAL STRAIN IN Z-DIRECTION IN COWL LIP; ELASTIC-VISCOPLASTIC ANALYSIS (ROBINSON’S MODEL). TIME. 0.75 sec; CROSS FLOW. 
GH 2 COOLANT. 


STRESS 

ksi 13.7 11.6 7.40 3.18 -1.04 -527 -9.49 -13.7 -17.9 



MPa 94.5 80.0 51.0 21.9 -72 -36.3 -655 -94.5 -123.4 



Z X 

FIGURE 6. - STRESS IN Z-DIRECTION IN COWL LIP; ELASTIC- VISCOPLASTIC ANALYSIS (ROBINSON’S MOOEL). TIME. 225 sec; CROSS FLOW. 
GH 2 COOLANT. 
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FIGURE 7. - CYCLIC THERMAL AND MECHANICAL LOADING USED TO 
GENERATE HYSTERESIS LOOPS BY THE MARC FINITE -ELEMENT 
CODE FOR FREED'S MODEL. 


FIGURE 8. - PREDICTED STABILIZED HYSTERESIS LOOPS USING 
FREED'S MODEL (ISOTHERMAL LOADING). EXPERIMENTAL DATA 
FROM FREED AND VERRILLI (1988). STRAIN RATE, O.OOI/SEC. 
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FIGURE 9. - IN-PHASE PREDICTED HYSTERESIS LOOPS FOR COPPER 
USING FREED'S MODEL. EXPERIMENTAL DATA FROM FREED AND 
VERRILLI (1988). TEMPERATURE, 200 Z 500 °C; STRAIN RATE, 
1.5x10~ 5 /sec. 


FIGURE 10. - OUT-OF-PHASE PREDICTED HYSTERESIS LOOPS FOR 
COPPER USING FREED'S MODEL. EXPERIMENTAL DATA FROM FREED 
AND VERRILLI (1988). TEMPERATURE, 200 Z 500 °C; STRAIN 
RATE, 1.5x10" 5 /SEC. 
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FIGURE 11. - CYLINDRICAL 
ROCKET NOZZLE FINITE- 
ELEMENT MODEL. ELEMENTS. 
35; NODES. 54. 
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FIGURE 12. - PROGRESSIVE SHAPE OF COMPONENT. PLOT 
MAGNIFICATION FACTOR. * 15. 
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FIGURE 13. - STRESS-STRAIN CURVES FOR SINGLE-CRYSTAL ALLOY 
(RENE-N4) . TEMPERATURE. 760 °C (1400 °F). 


FIGURE 14. - CREEP CURVES FOR SINGLE-CRYSTAL ALLOY (RENE-N4). 
DIRECTION (1 0 0); TEMPERATURE. 760 °C (1400 °F). 
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FIGURE 15. - RELAXATION CURVES FOR SINGLE-CRYSTAL ALLOY 
(RENE-N4). TEMPERATURE, 760 °C (1400 °F), 



STRAIN 

FIGURE 16. - HYSTERESIS LOOPS FOR SINGLE-CRYSTAL ALLOY 
(RENE-N4). DIRECTION (1 10); TEMPERATURE, 760 °C 
(1400 °F). 
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